ADA045535 


ryi^FTTI/M  IVtXX  FT  IE  T~rxx 


Martin  Marietta 
Laboratories, 


^{jioAir/  LA'iiai  ovHR  bluni’-body 

OF  PB'OIin’I'jN  AT  INCIDENCE  A'TJ  NAGNUS 
F0RC2S 


{[.  C.  l*S.Tlg 


/vii'SiV.st  19'^’r 


O#' 


Rfic  larch  siionsorofi'.  by  the  IJ.  S.  Am>- 
Re^<;arch  Office  urder  Cont'act  !lM-CC4-7.‘> 
C'OOll  end  by  the  U.  S.  M:'  Forci*  C’t'fice 
of  Scie'.itific  Fjesearrfi  ^.’PSC)  > 

Contract  F49d20-76-C-000d . 


I DISTRIBin’lON  STATHM^iNT  A 

! Appioviid  Ici  pi'.tiio  isjk'ciso; 

! ihiihuCcd 


BOUNDARY  LAYER  OVER  SPINNING  BLUNT- 
BODY  OF  REVOLOTION  AT  INCIDENCE  AND  MAGNUS  FORTES 


K.  C.  Wang 


NML  TR  77-SOc 


August  1977 


Research  Sponsored  by  the  U.  S.  Array  Research 
Office  under  Contract  DAHC04-73-C-0011  and  by 
the  U.  S.  Air  Force  of  Scientific  Research  (AFSC) 
under  Contract  F49620-76-C-0004. 


Approved  for  Public  Release  - Distribution  Unlimited 


Martin  Marietta  Corporation 
Martin  Marietta  Laboratories 
1450  South  Rolling  Road 
Baltimore,  Maryland  21227 


SKCUNITy  classification  of  this  page  r***>«>  Omtm  Eniarvo 


REPORT  DOCUMENTATION  PAGE 


ITLC  (md  Subili/*; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT’S  catalog  NUMBER 


^xmdary  ^yer  Jgpinning 
Revolution  at  incidence  and^ 


MtNG  OACANIZATIOM  NAME  AND  AOOACSS 

Martin  Marietta  Laboratories 

Martin  Marietta  Corporation  ^ 

1450  South  Rolling  Road,  Baltimore,  Maryland  21227 


NML  TR  77-SOc 


8.  contract  or  grant  NUMBERT*; 


DAHC04-73-C-^l 

F496;0-76-C-Bf00^ 


n.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 


13.  NUMBER  OF  PACES 

83 


ITORING  agency  NAME  S AODRESSCJ/  di/<*rant  front  Coniratlini  Otffeo)  | IS.  SECURITY  CLASS,  (of  t/i/t  nportj 

Unclassified 


fmcy 


ISa.  declassification/oowngraoing 

SCHEDULE 


IS.  distribution  statement  (of  thin  Rnport) 


.^proved  for  public  release;  distribution  unlimited. 


t7.  DISTRIBUTION  STATEMENT  (oi  In  BlocA  2C,  it  dUf^rmit  from  Kopcrt 


IS.  KEY  WORDS  (Contlrtum  on  rororoo  oldo  it  nocoooofr  md  idontiiy  by  block  ntmboe) 


Three-dimensional  boundary  layer 
Boundary  layer  over  spinning-body  at  incidence 
Magnus  forces 
Numerical  solutions 


abstract  (Contlnao  m rooofoo  aido  it  nocooomr  md  Idmtity  by  bfoeir  nuoibor) 

An  incompressible  laminar  flow  over  a spinning  blunt -body  at  incidence  is  in- 
vestigated. The  approach  follows  strictly  the  three-dimensional  boundary  layei 
theory,  and  the  lack  of  initial  profiles  is  readily  resolved.  The  rule  of  the 
dependence  zone  is  satisfied  with  the  Krause  scheme  and  con^lete  numerical 
solutions  are  obtained  for  an  ellipsoid  of  revolution  at  6*  , incidences  and  two 
different  spin  rates.  Spinning  causes  asymmetry  vdiich,  in  tmn,  introduces 
the  Magnus  force.  The  asyiimetry  is  most  pronounced  in  crossflow,  but  is 


Farrm 
I JAN  TJ 


edition  of  I NOV  «■  IS  OBSOLETE 


Unclassified  ^ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  Onto  Kntmtnd) 


SKCUNITY  CUASSIPICATION  op  this  PAOC(l*>l«n  OMa  Cnttn^ 


^noticeable  in  the  skin-friction  and  displacement  thickness  of  the  meridional 
flow.  A variety  of  crossflow  profiles  are  determined  as  are  the  streaml^e 
patterns  in  the  cross-  and  meridional -planes  which  are  espci^ly  i^e^  in 
visiaalizing  the  flow  structure. Detailed  distribution  of  skin  frictions,  dis- 
placement thicknesses,  and  centrifiigal  pressxnre  are  presented.  A nw 
derivation  of  the  total  displacement  thickness  is  given  and  a negative  dis- 
placement thickness  is  found  physically  meaningful  The  Magnus  forces  due 
to  the  crossflow  skin  friction  and  the  centrifugal  pressure  are  determined; 
these  two  forces  partly  compensate  each  other.  At  lower  spin  rage,  the 
frictional  force  is  larger,  resulting  in  a positive  Magnus  force.  At  high  spin 
rate,  the  opposite  is  obtained.  At  high  incidence  (30°),  the  leeside  separated 
region  associated  with  an  open  separation  is  found  not  amenable  to  a classical 
boundary  layer  treatment.  The  present  boundary  layer  calculations  could  be 
carried  out,  in  the  longitudinal  direction,  only  up  to  the  beginning  of  an 
open  separation.  Since  an  open  separation  moves  forward  with  increasing  inci- 
dence, the  calculable  area,  therefore,  decreases. 


wns  tv.  o' 


SCCUKITY  classification  of  this  NAOErWhan  Dtm  SnIcrMi) 


CONTENTS 


» 


ABSTRACT 

^ 1.  INTRODUCTION 

2.  FORMULATION 

^ 2.1  Geometry  and  Coordinates 

2.2  Boundary  Layer  Equations 

2.3  Boundary  Layer  Thicknesses 

^ 2.4  Magnus  Forces 

2.4  a.  Sectional  Friction  Side-Force 

2.4  b.  Sectional  Centrifugal  Side-Force 

I 

3.  MEfflOD  OF  SOLUTION 

3.1  Boundary  Layer  Method 

3.2  Finite  Difference  Schemes 

3.3  Computation  Details 

3 . 4 Starting  Method 

I 3.4  a.  Initial  Profiles  for 

Meridional  Marching 

3.4  b.  Initial  Profiles  for  Circumferential 

Marching 

^ 3.5  Comparison  With  Other  Methods 


CONTENTS 


4.  RESULTS  AND  DISCUSSIONS  28 

4.1  Inviscid  Input  28 

4.2  Boundary  Layer  Profiles  28 

4.3  Cross -plane  Flow  Pattern  32 

4.4  Meridional -Plane  Flow  Pattern  39 


4.5  Meridional  (Or  Longitudinal)  Skin 
Friction, 

4.6  Circumferential  (or  Cross-flow)  Skin 


Friction,  C£g  42 

4.7  Centrifugal  Pressure,  p^^  44 

4.8  Displacement  Thicknesses  46 

4.8  a.  v-Velocity  Displacement 

Thickness,  A*  46 

* 

4.8  b.  Total  Displacement  Thickness,  A 48 

4.9  Magnus  Forc§  50 

4.10  The  High- Incidence  Problem  52 

5.  CONaUSIONS  58 

ACKN0WLEDGB4EOTS  60 

6.  REFERENCES  61 

APPENDIX  A - Derivation  of  the  Displacement 

Thickness  Equation  63 

APPENDIX  B - Additional  Figures  67 


ABSTRACT 


An  incompressible  laminar  flow  over  a spinning  blunt -body  at 
incidence  is  investigated.  The  approach  follows  strictly  the  three- 
dimensional  boundary  layer  theory,  and  the  lack  of  initial  profiles 
is  readily  resolved.  The  rule  of  the  dependence  zone  is  satisfied  with 
the  Krause  schane  and  complete  numerical  solutions  are  obtained  for  an 
ellipsoid  of  revolution  at  6®  incidences  and  two  different  spin  rates. 
Spinning  causes  asyranetry  which,  in  turn,  introduces  the  Manus  force. 

The  asyirmetry  is  most  pronomced  in  crossflow,  but  is  also  noticeable  in 
the  skin-friction  and  displacement  thickness  of  the  meridional  flow.  A 
•variety  of  crossflow  profiles  are  determined  as  are  the  streamline 
patterns  in  the  cross-  and  meridional -planes  which  are  especially  useful 
in  -visualizing  the  flow  structure.  Detailed  distribution  of  skin  frictions^ 
displacement  thicknesses,  and  centrifugal  pressure  are  presented.  A neiv 
derivation  of  the  total  displacement  thickness  is  given  and  a negative 
displacement  thickness  is  found  physically  meaningful.  The  Magnus  forces 
due  to  the  crossflow  skin  friction  and  the  centrifugal  pressure  are 
determined;  these  two  forces  partly  compensate  each  other.  At  lower  spin 
rate,  the  frictional  force  is  larger,  resulting  in  a positive  Magnus  force. 
At  high  spin  rate,  the  opposite  is  obtained.  At  high  incidence  (30®), 
the  leeside  separated  region  associated  with  an  open  separation  is  found 
not  amenable  to  a classical  boundary  layer  treatment.  The  present 
boundary  layer  calculations  could  be  carried  out,  in  the  longitudinal 
direction,  only  up  to  the  beginning  of  an  open  separation.  Since  an  open 
separation  moves  forward  with  increasing  incidence,  the  calculable  area, 
therefore,  decreases. 


1 


INTRODUCTION 


The  three-dimensional  laminar  boundary  layer  over  an  inclined 
body  of  revolution  spinning  about  its  awi  major  axis  is  considered. 

Aside  from  interest  in  its  own  right  as  a moving-wall  boundary 

layer  flow  in  three  dimensions,  this  kind  of  problem  has  long  attracted 

attention  in  connection  with  the  "Magnus  effect"  for  a flying  projectile. 

The  idea  of  Magnus  effect  is  usually  introduced  in  the  literature  as 


the  lift  force,  prV  , on  a rotating  circular  cylinder CFig-  la). 


It  is  presented  within  the  context  of  a potential  flow  in  terms  of  a 


pry. 


Ca) 


Fig.  1.  Magnus  Force. 


uniform  velocity  and  an  assumed  circulation  F,  although  the  circulation 
was  noted  to  be  due  to  the  viscous  effect.  However,  the  precise  connection 
betwen  the  circulation  and  the  viscous  flow  (or  boundary  layer  flow) 
has  always  been  difficult  to  make  very  clear. 

For  an  inclined  finite  long  body  (Fig.  lb),  the  Magnus  force  is 
defined  to  be  the  side-force  perpendicular  to  the  lift  plane  (containing 
the  lift  vector  and  the  free  stream  velocity).  The  mcmient  due  to  such 


a side -force  is  known  as  the  Magnus  moment.  The  Magnus  force  and 
moment  may  cause  a flying  object  to  depart  from  its  pre-set  tra- 
jectory. The  cross-sectional  view  of  Fig.  lb  resembles  the  two- 
dimensional  picture  of  Fig.  la,  but  this  way  of  looking  at  the  problem 
doesn't  really  help  much.  In  a three-dimensional  problem,  the  circu- 
lation is  not  always  a very  convenient  concept,  and  the  Magnus  forces 
could  be  better  analyzed  in  terms  of  surface-distributed  forces.  The 
question  about  what  surface  force  would  contribute  to  the  Magnus 
effects  and  their  calcxilability  turns  out  to  be  a very  challenging  one. 
Our  search  to  understand  the  mechanism  and  to  determine  quantitatively 
the  Magnu'  ;s  continues  to  the  present. 

Mr  was  the  first  to  consider  the  three-dimensional  boundary 

layer  over  a cylindrical  portion  of  a body  of  revolution  at  small  inci- 
dence and  spin  rate  by  a small  perturbation  method  and  to  identify  the 
asymmetry  of  boundary  displacement  thickness  as  a source  of  Magnus  force. 
A similar  approach  was  independently  carried  out  by  Kelley . Sedney^^^ 
further  extended  Martin's  analysis  to  the  case  of  a supersonic,  sharp 
cone.  These  investigations  correctly  point  out  the  right  direction  for 
latter  work,  even  through  the  results  were  subject  to  the  limitations 
of  small  perturbations. 

Interest  in  this  problem  has  resurfaced  in  recent  years  and  the 
method  of  small  perturbations  has  been  replaced  by  the  finite- 


difference  method,  thus  removing  the  restrictions  of  small  incidence  and 
small  spin.  Dwyer^^^  led  the  way,  followed  by  Clark^^^ , Vaughn  and 
Reis^^\  Watkins Lin  and  Rubin^^^  and  Dwyer  and  Sanders 
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Dwyer calculated  the  three-dimensional  boundary  layer  over  a 
hypersonic  sharp  cone,  but  encountered  difficulty  in  obtaining  a 
unique  solution.  This  led  him  to  conclude  that  the  boundary  layer 
equations  are  invalid  for  the  spinning-body  problem,  and  that  the 
lateral  diffusion  term  should  be  retained  in  the  governing  equations. 

Dwyer's  contention  was  later  shared  by  Clark and  by  Lin  and 

Rubin^^^ . Clark  adopted  full  Navier-Stokes  equations  for  the  case  of 

a spinning  sphere  and  an  ellipsoid  of  revolution  at  incidence,  but  lack 

of  sufficient  resolution  apparently  presents  a difficult  problem  and 

fg") 

very  sketchy  results  for  the  ellipsoid  were  given.  Lin  and  Rubin 
employed  a parabolicized  version  of  the  Navier-Stokes  equations  similar  to 
that  suggested  by  Dwyer ; they  obtained  a coirplete  solution  for  the 
spinning  cone  problan  and  singled  out  the  circumferential  skin  friction 
as  a major  contribution  to  the  Magnus  force. 

Watkins^ disputed  Dwyer’s  conclusions  and  attenq)ted  to  calculate  the 
boundary  layer  flow  over  a spinning  ogive -cylinder  at  incidence.  Watkins' 
insistence  on  following  the  boiindary  layer  approach  is  correct,  but  there 
are  apparently  other  difficulties  which  make  his  work  incon^Jlete  with 
scanty  results.  Voughn  and  Reis^^^  developed  an  approximate  analytical 
solution  for  spinning  projectiles.  Their  method  should  be  particularly 
attractive  for  engineering  design  purposes.  These  authors  included  the 
centrifugal  pressure  as  another  important  contribution  to  the  Nfeignus  force. 
Dwyer later  reconsidered  the  spinning  cone  problem  in  collaboration 
with  Sanders.  This  time,  they  returned  to  the  boundary  layer  approach,  and 
unique  solutions  were  obtained.  In  particular,  the  Maenus  force  due  to  the 


I 


[ 

J 
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aisynnetric  displacement  thickness  was  calculated  by  a separate 
inviscid  program  in  contrast  to  the  slender  body  theory  used  by 
previous  authors However  their  modified  difference  scheme 
still  does  not  completely  satisfy  the  rule  of  dependence  zone.  This 
will  be  further  discussed  later  in  Sect.  3.5. 

Spinning  causes  reversal  of  circumferential  flow  and  whether 
calculation  of  such  reversed  flow  would  violate  the  initial -value  con- 
cept for  the  boundary  layer  equations  has  been  an  open  question.  Lack 
of  a strong  theoretical  basis  contributes  to  the  confusion  about  the 
adequacy  of  boundary  layer  theory  and  the  necessity  of  using  parabola cized 
or  full  Navier-Stokes  equations.  Of  course,  higher  approximation  offers, 
in  principle,  a more  complete  flow  description  than  the  boundary  layer 
theory  can  pixavide.  However  for  specific  purposes,  this  may  not  always 
be  necessary  nor  desirable.  More  elaborate  fomnulation  requires  longer 
computing  time,  and  may  lead  to  insufficient  resolution,  making  results 
unreliable  and  inconclusive. 

In  related  work^^^^,  a similar  reversed  circumferential -flow  was 
calculated  by  this  author  following  strictly  the  three-dimensional 
boundary  layer  formulation.  This  experience  convinces  us  that  the  spinning 
problem  can  be  similarly  treated. 

In  this  work,  we  present  our  detailed  methods  and  results.  The 
specific  problem  investigated  is  concerned  with  an  incompressible  flow 
over  a spinning  ellipsoid  of  revolution  at  incidence.  The  same  problem 
without  spinning  was  considered  earlier . Such  a blunt-body 
problem  provides  interesting  contrasts  to  the  supersonic  sharp  cone  case 
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previously  studied  because  of  the  presence  of  meridional  (longitudinal) 
pressure  gradients,  favorable  and  adverse.  Lack  of  pressure  gradients 
would  mean  less  boundary- layer  variation.  Section  2 describes  the  formu- 
lation, including  a new  rederivation  of  the  total  displacement  thickness 
considered  by  Moore  and  Lighthill  . Section  3 explains  the  method 

of  solution,  the  rationale  of  reversed-flow  calculation,  and  the  starting 
procedure.  Lack  of  initial  profiles  was  considered  as  an  insurmountable 
difficulty  for  the  boundary  layer  approach,  but  actually  this  question 
can  be  readily  resolved  as  demonstrated  here.  Section  4 presents  thf  re- 
sults and  discussions.  Complete  distributions  of  various  boundary  1 «.yer 
diaracteristics  on  the  spinning  body  are  presented.  In  contrast,  results 
of  the  distribution  of  skin  frictions,  for  exaii5)le,  for  the  pointed  bodies 
previously  studied  are  either  completely  absent or  only  given  at  a 
single  crossplane Section  5 gives  the  conclusions,  j^pendix  A 
details  the  derivation  of  the  equations  for  total  displaceir-ent  thckness 
in  three-dimensional  boundary  layers.  Appendix  B provides  more  calculated 
results  than  those  included  in  the  main  text. 

Highlights  of  this  work  were  briefly  reported^^^^  before  with  re- 
stricted distribution;  the  preparation  of  the  present  report  has  been 
interrupted  and  delayed.  This  work  was  presented  at  the  14th  International 
Congress  of  Theoretical  and  Applied  Mechanics  at  Delft,  September  1976. 

More  recent  calculations  have  extended  the  same  problan  to  the  high  inci- 


dence case. 
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2.  FORMULATION 
2.1  Geometry  and  Coordinates 

fl21 

As  in  our  previous  study'-  , the  specific  example  of  an  ellipsoid 
of  revolution  (Fig.  2)  vdth  minor/major  axes  ratio  of  1/4  was  chosen 
for  concrete  results.  A non-rotating  body  coordinate  system,  p,  9,  z is 


Fig.  2.  Geometry  and  Coordinates. 

used.  As  depicted  in  Fig.  2,  a is  the  incidence  angle,  p and  0 are 
two  surface  coordinates  along  the  meridional  and  circumferential 
directions,  respectively,  and  z denotes  the  noimal  coordinate.  The 
corresponding  metric  coefficients  are 

h^^  = [(l-eV)/(l-P^)]^, 

hQ  = [(l-e2)Cl-p2)] 
vAiere  e is  the  eccentricity  defined  by 

e » [l-(b/a)^]**, 

a and  b being  the  semi -major  and  semi -minor  axes.  h.  is  just  the  cross 
sectional  radius  at  a fixed  p- station.  We  shall  denote  u,  v,  w as  the 


velocities  along  n,  0,  z directions;  u and  v are  non-dimensionalized 

u 

with  the  flow  velocity  at  infinity  V^;  w with  V^yR  vrtiere  RC“V^ 

is  the  Reynolds  number;  z is  non-dimensionalized  with  a/R  ; the 
2 

pressure,  with  pV^  , \diere  p is  the  density.  The  ellipsoid  spins 
about  its  major  axis  with  an  angular  velocity  <*>  non-dimensionalized 
with  V^a.  The  direction  of  spin  is  clockwise  viewed  from  the  rear. 
All  the  variables  are  in  the  physical,  primitive  form;  i.e.,  no  trans- 
formed variable  is  used. 

2.2.  Boundary  Layer  Equations 

A strict  boimdary  layer  approach  is  followed  here.  The  govern- 
ing equations  in  the  coordinate  system  of  y,  9,  z are 


..  3u 

3u 

+ w M - v^- 

^ . 

, CD 

TTly 

3z^ 

3v 

3v 

3v 

^ -UT  + 11V 

3h. 

, ■ P s 

+ ^ . 

(2) 

h hQ3y 
y 0 

3z 

3u 

3v 

3z 

- 0, 

(3) 

with  the  boundary  conditions 

u*W"0, 

V = uhg  at 

z*0. 

(4a) 

u»u,  v»V 

at  z-*«. 

(4b) 

Ik 


v*ere  U and  V are  the  velocities  at  the  boundary  layer  edge.  In  the 
present  case  of  a prolate  spheroid,  U and  V are  known  from  an  exact 
inviscid  flow  solution , 

^ * 71  ri!.  U'  Ccosa)  (1-y^)**  + 0/a)  (1+kJ  (sina)ycosel , (Sa) 

Cl-e  vj  *"  ^ -I 

V = sina  sin0,  (5b) 

where  k^  and  k^  are  known  as  the  axial  and  cross -coefficients  of 
virtual  mass  and  are  defined  by 

‘'a  • i^]  • 


1 

1+21C- 


Knowing  U and  V,  the  pressure  is  given  by 


l-(u2.v^). 


The  two  components  of  skin  friction,  meridional  c^  and  circumferential 

2 

c^g,  non-dimensionalized  with  pV^  are  defined  by 


- — 

ft  ^ 

(m 

\ 3Z/Z-K)  » 

C8a) 

n . 1 i 

f 3v\ 

Jz)z-*0  * 

(8b) 
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2.3.  Boundary  Layer  Thicknesses 


The  two  cori5X3nents  of  displacanent  thickness,  A*  and  A ^ , 


I ised  respectively  on  the  u-  and  v-velocity  profiles  are  defined  by 

00 

and 

1 r 

= i/CV-v)d2,  C9b) 


where  q is  the  total  inviscid  velocity, 

q = (U^  + . 


The  four  components  of  momentum  thickness,  based  on  u-  and  v- 
velocity  profiles  and  their  cross  products,  are  defined  by 


*■  7) 


GJ-u)udz, 

(10a) 

) 

CV-v)vdz, 

(10b) 

> 00 

GJ-u)vdz  , 

1 

(10c) 

00 

(V-v)udz  . 

(lOd) 
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Notice  that  these  boundary  layer  thicknesses  are  deliberately  not 
written  in  non-dimensional  form.  This  was  done  to  avoid  the  situation 
vdien  the  inviscid  velocity,  V,  becomes  zero  at  0 = 0°  and  180°.  The 
vanishing  of  V does  not  imply  the  vanishing  of  v-profile  in  a moving - 
wall  problem. 

The  displacement  thickness  of  a three-dimensional  boundary  layer 
A*  is  defined  in  terms  of  A*  and  A*  by  the  equation: 

^ [hg  OJA*-  qA^)]  + Ig-  [hy  (VA*  - qAg*)]  = 0.  (11) 

This  was  first  derived  by  Moore Later  Lighthill gave  four 
alternative  methods  of  derivation.  Moore's  original  method  was  included 
and  termed  by  Lighthill  as  the  "velocity  comparison"  method.  Lighthill's 
derivation  relies  on  the  streamline  coordinates  formed  by  the  equi- 
potential  lines  and  the  inviscid  streamlines,  vdiereas  Moore's  derivation 
is  valid  in  general  orthogonal  coordinates.  In  Appendix  A,  a new  alterna- 
tive derivation  is  presented.  Our  derivation  is  patterned  after  that  of 
the  ordinary  mass  continuity  equations;  the  reasoning  is  somewhat  akin 
to  Lighthill's  "flow  reduction"  method,  but  differs  in  the  formalism. 

Also  it  does  not  depend  on  the  particular  coordinate  systan.  We  shall 

* * 

refer  to  A as  the  total  displacement  thickness,  A^  as  the  u-velocity 

(or  meridional  flow)  displacement  thickness,  and  Ag  as  the  v-velocity  (or 
crossflow)  displacement  thickness. 

2.4  Magnus  Force 

The  Magnus  force  results  from  the  flow  asymmetry  \diich,  in  turn,  is 
caused  by  the  spinning  motion  through  the  action  of  fluid  viscosity. 
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Symmetry  refers  here  to  the  symmetry  plane  of  the  inviscid  flow  or  the 
plane  of  9 * 0°  and  180°  of  Fig.  2.  As  mentioned  in  the  Introduction, 
three  major  contributions  to  the  Nfagnus  force  have  been  identified: 

Cl)  the  circumferential  skin  fiction,  (2)  centrifugal  pressure  and  (3) 
pressure  due  to  the  asynmetrical  displacement  thickness.  Only  the  first 
two  are  determined  here;  the  third  one  requires  a separate  inviscid  ccan- 
puting  program  and  hence  has  not  been  calculated.  It  has  been  known  that 
the  asynmetric  wake  vortices  shed  from  an  inclined  body  of  revolution 
even  for  a non-spinning  case  also  generate  side-forces.  However  deter- 
mination of  these  wake  vortices  requires  a conplete  Navier-Stokes  solution 
and  hence  is  beyond  the  scope  of  the  present  work. 

2.4a.  Sectional  Friction  Side-Force 


The  sectional  side-force  due  to  the  circumferential  skin  friction  is 
given  by 


Fig.  3a.  Frictional  side-force. 

2 

where  F^^  is  nondimens ionali zed  with  pV^  , positive  to  the  right.  The 
subscript  sf  stands  for  skin  friction.  A corresponding  force  coefficient, 

Csf,  is  defined  by 


>/R 


>/R  F, 


Sf 


>/R  F, 


sf 


sf 


area 


irCb/a) 


(12b) 


i.e.,  the  reference  area  is  taken  to  be  the  mid-section  = 0) 
of  the  ellipsoid. 

2.4.b.  Sectional  Centrifugal  Side-Force 

The  sectional  side-force  due  to  centrifugal  effects  is  ob- 
tained from  the  centrifugal  pressure  p^^  vdiere  the  sxobscript  ct  stands 
for  "centrifugal."  After  v is  known  p^^  is  evaluated  from, 

H ■ ^ 


Replacing 


r - h.  + -—  Cos  ((» a h.  , 
® n/R  ® 


^ ^ 3p  3z  _ ^ 

5r  * "Sz  "Sr  ~ dz  cos  <{) 


Pv  + p.,^/Vr  , 

body  edge  ct' 


Fig.  3b.  Centrifugal  pressure. 


\diere 


C13a) 


Integrating  this,  we  obtain 


Fig.  3c.  Centrifugal  side-force. 


where  F^^  again  is  taken  positive  to  the  right.  Substituting 

^ ^ct  “ "/  [■  r* ^ ^ ’ 

.-if  ^^9  \ 

■ sin  I (TTrr  J • 


(13b) 
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^edge  is  the  inviscid  surface  pressure.  It  is  unaffected  by  the 
spinning  motion  and  is  synmetrical  with  respect  to  the  plane  of  9»0° 
and  180°.  Therefore  integration  around  the  body  vanishes.  The 
corresponding  centrifugal  force  coefficient  is 
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3.  MEIHDD  OF  SOLUTION 
3.1  Boundary  Layer  Method 

The  boundary  layer  method  is  followed  through  in  this  work. 

This  choice  comes  naturally  because  of  our  previous  experience  in  a 
non-spinning  case  for  the  same  body  where  similar  reversal  of  the 
circumferential  flow  also  occurs.  The  calculability  of  such  reversed 
flow  was  rationalized  thereby  using  the  rule  of  the  dependence  zone 
pertaining  to  three -dimensonal  boundary  layer  equations.  What  are 
really  involved  can  be  seen  from  Fig.  4,  which  depicts  the  dependence 


Fig.  4.  The  zone  of  dependence,  (a)  coordinates;  Cb)  nonspinning, 
no  reversal;  (c)  nonspinning,  reversal;  (d)  spinning,  no 
reversal;  (e)  spinning,  reversal. 

zone  (shaded  area)  in  the  p-0  plane  (i.e.  the  body  surface)  for  a 
point  P.  Fig.  4a  indicates  the  coordinates,  assvmdng  that  the  body 
spins  along  the  positive  9 -direction.  Fig.  4b  indicates  that  when 
there  is  no  flow  reversal,  the  dependence  zone  of  point  P falls  into 
the  third  quadrant.  When  the  non-spinning  circumferential  flow  reverses, 
the  dependence  zone  situates  across  the  y-coordinate  line  (Fig.  4c). 
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Spinning  merely  makes  the  lovfer  edge  of  the  dependence  zone  almost 
coincide  vdth  the  9 -coordinate  lines  (Fig.  4d.e).  As  long  as  the 
rule  of  the  dependence  zone  is  satisfied  by  the  difference  scheme 
used,  there  are  no  logical  differences  among  these  cases.  Reversal 
of  one  component  of  the  flow  is  not  the  same  as  reversal  of  the  total 
flow,  and  hence  does  not  automatically  lead  to  violation  of  initial -value 
idea.  In  fact,  the  flow  direction  (parallel  to  the  body  surface) 
varies  across  the  three-dimensional  boundary  layer  at  a given 
point  on  the  body,  thus  sweeping  out  a certain  solid  angle. 

A distinctive,  single  resultant  flow  direction  does  not  exist  in  this 
case.  The  dependence  rule,  vdiich  requires  the  crarputation  mesh  to 
enclose  the  solid  angle,  ensures  precisely  the  satisfaction  of  the 
initial -value  concept. 

3.2.  Finite-Difference  Scheme. 

The  finite  difference  schemes  used  in  this  work  are  those  employed 
previously  for  the  non-spinning  case;  they  are  implicit  schemes  of  the 
Crank-Nicolson  type.  Corresponding  to  the  dependence  zones  of  Fig.  4, 
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Three  different  meshes  are  used  (Fig.  5).  For  convenience,  these  meshes  will 
be  referred  to  as  schemes  1,  2 and  3.  The  shaded  angle  denotes  each  scheme's 
dependence  angle.  Scheme  1 is  used  for  area  between  0=0®  and  180®  where 
no  negative  circumferential  flow  occurs;  scheme  2 is  used  for  generating 
initial  profiles  along  the  plane  9 = 0°  or  360°,  which  will  be  further 
discussed  later  (Sect . 3 . 4b) . Scheme  3 is  used  for  the  area  with  both 
positive  and  negative  circumferential  flows.  The  data  at  the  solid 
points  are  used  to  calculate  solutions  at  the  ai^jty  circles  in  a single 
calculation  step.  Typical  approximations  for  scheme  1 are 


3u^  2Am  ^A+l,m+l,n'^  ^Jl,m+l,n  ^ ^ll+l,m,n  " ^ ^£,m,nj  ’ (15a) 


30  2A0  ^A+l,m+l,n  ^ ^A+l,m,n  * ^ ^A,m+l,n  ' ^ ^A,m,nJ  ’ 


(15b) 


For  schene  2; 


■ ISy  [ ^ h+l,m,n  » (16a) 

2a9  ^Z,m+l,n  ' ^ H,m-l,nj  ’ (16b) 

~5i^“  “?AF  [f  ^A+l,m,n+l  ^A+l,m,n-l  * ^ ^i-l,m,n+l  ' ^ ^i-l,m,n-l]  » C16c: 


"iz^  2(Az)^  *■  ^^4+l,jn,n  ’ ^ ^^Z-l,ni,n]' 


(16d) 


For  scheme  3; 


Sy  ^ ~ Ay  [ ^ ^£+l,m,n  ' ^ ^i,m,n]  ' 


(17a) 


lU. 

36 


Sr  ^i+l,m,n  ■ ^ ^Z+l,m-l,n  ^ ^Jl,m+l,n  " ^ ^Jl,in,n]  ’ 


(17b) 


^ ^i+l,m,n+l  ■ ^ ^A+l,m,n-l'^  ^ ^il,m,n+l  ' ^ ^i,m,n-l]’ 


W,n-C^^C)),,,J. 


(17d) 


For  the  calculation  of  w from  the  continuity  equation,  the  derivative 
3w/3z  is  approximated  at  a single  central  station  in  these  respective 
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Fig.  6.  Difference  meshes  for  continuity  equation. 


schemes  (Fig*  6).  As  discussed  in  Ref.  12,  this  derivative  becomes 


dz 

('^)l+35,m+Js,n 

(Scheme  1) 

(18a) 

^ =1 
3z  ' 

^'^Z,m,n  Az  , 

(Scheme  2) 

(18b) 

=> 

('^i+Js,m,n 

'^il+J5,m,n-l)/^  ' 

(Scheme  3) 

(18c) 

At  first,  this  derivative  was  approximated  by  averaging  over  comer 
stations  and  was  found  to  exhibit  fliictuation,  causing  eventual  breakdown 
of  calculations.  Such  difficulty  results  frcan  the  difference  approximation 
for  the  continuity  equation  and  hence  falls  beyond  the  scope  of  usual  nimieri 
rai  Stability  analysis  applied  to  the  linearized  version  of  the  momentum 
equations.  Linearization  decouples  or  separates  the  instability  problem 
between  the  momentum  equations  and  the  continuity  equation.  At  the  present 
time,  there  is  no  known  method  for  conducting  a numerical  stability  analysis 
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for  a system  of  nonlinear  coupled  equations,  and  hence  the  difficulty 
in  calculating  w just  mentioned  has  not  been  clearly  understood. 

Scheme  1 was  first  used  by  Dwyer scheme  2 by  Hall^^^^ 
and  scheme  3 by  Krause They  were  previously  employed  by  the 
present  author  for  the  corresponding  nonspinning  cases 

3.3  Conq3utation  Details 

As  in  the  nonspinning  cases,  our  con^jutation  proceeds  first  along 
the  circumferential  direction  from  9 = 0 to  it  and  to  2it  , and  then  aft 


Fig.  7.  Computation  march. 

along  the  meridional  direction.  Scheme  1 was  used  between  9=0  and^f 
where  there  is  no  circumferential -flow  reversal.  Scheme  3 was  used  be- 
tween 9 * IT  and  2ir  where  circumferential  reversed  flow  always  prevails. 

It  should  be  noted  that  reversed  circumferential -flow  also  occurs  at  9 
less  than  180®  (say  120®  -v  180®)  for  the  nonspinning  case,  but  such 
reversed  flow  is  annulled  by  opposite  flow  due  to  spinning.  This  is  tiue 
at  least  for  the  case  calculated  here.  Schane  1 has  the  advantage  of 
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! 1 

I siji5)licil:y  and  is  available  from  our  previous  development.  Other - 

i wise  scheme  3 can  alternatively  be  used  for  the  vdiole  calculation. 

The  conqniter  program  has  been  built  to  allow  different  step-sizes 
in  the  normal  direction  z , but  only  two  different  steps  were  used  in 
the  present  calculations.  Smaller  step  was  maintained  near  the  body 
surface  to  ensure  more  accurate  evaluation  of  skin  frictions,  although 
the  results  show  that  skin  frictions  are  not  really  much  affected  by 
the  change  of  Az.  On  the  average,  50  steps  were  taken  in  the  z -direction. 

In  order  to  save  coi^uting  time^  experiments  of  using  different  step- 
sizes  along  the  y-  and  9-direction  were  made.  Larger  steps  can  be 
employed,  depending  on  how  severe  the  relevant  adverse  pressure  gradient 
is.  The  case  of  cu  = 1.0  at  6®  incidence  was  calculated  with  Ay  = 0.05 
and  0 = 10®,  and  the  vdiole  calculation  took  less  than  one  hour  of  CPU 
time  on  a UNTVAC  1108  con^juter. 

3.4  Starting  Ifethod 

To  start  the  calculation,  initial  profiles  are  needed  along  both 
surface  coordinate  directions.  For  a nonspinning  inclined  body  of 
revolution,  a separate  boundary-layer  solution  along  the  windward  symmetry 
plane  has  often  been  used  to  provide  the  required  initial  profiles  for 
conpjtation  along  the  circumferential  direction.  For  the  corresponding 
spinning  case,  lack  of  such  a symmetry  plane  and  hence  its  associated 
boundary- layer  solution  has  caused  a major  difficiUty. 

Initial  profiles  for  boundary  layer  calculation  generally  do  not 


have  to  be  exact.  The  boundary- layer  problem  is  a mixed  initial-  and  boundary- 
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value  problem.  The  boundary-value  part  of  the  problem  provides  some 
capability  for  self -correcting  reasonably-approximated  initial  profiles. 

The  adjective  "reasonably  approximated"  must  be  emphasized.  If  the 
initial  profiles  are  arbitrarily  picked,  the  calculation  may  be  doomed, 
particularly  for  three-dimensional  problems,  in  comparison  to  two- 
dimensional  ones. 

Lack  of  synmetry-plane  solutions  led  Dwyer to  resort  to  a 

"guess"  method  with  iteration,  which  was  considered  to  be  semi -sat is factory, 

except  that  only  some  oscillations  in  the  results  were  found  near  the  plane 
f221 

9=0.  Watkins^  ^ carried  out  a small  perturbation  solution  about  the 

nonspinning  synmetry-plane  solution  that  is  valid  for  moderate  spin 

rates.  This  was  used  to  develop  initial  profiles  in  his  later  work^®^. 

For  large  spin  rates,  the  same  solution  was  suggested,  coupled  with  cyclic 

iterations.  Since  only  sketchy  results  were  presented,  it  is  difficult  to 

judge  how  this  starting  method  worked  out. 

(9  10  231 

Other  authors^  ’ ’ ''  sought  to  develop  difference  schemes  that  do 

not  require  such  circumferential  initial  profiles.  These  alternative 
methods  have  objectionable  side -effects.  The  starting  methods  used  here 
will  demonstrate  that  such  initial-value  questions  can  be  satisfactorily 
resolved.  Comparisons  with  the  above-mentioned  methods  will  be  made  later. 

3. 4. a.  Initial  Profiles  for  Meridional  Marching 

Initial  profiles  along  a fixed  crossplane  near  the  vertex  of  the 
body  are  needed  for  meridional  march  of  conputation.  Near  the  vertex, 

spinning  effects  have  not  yet  materially  changed  the  corresponding  non- 
spinning profiles  except  for  the  v-velocity  profile  near  the  body  surface. 
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I Based  on  this  observation,  initial  u-profiles  were  taken  to  be  the 

1 same  as  those  of  the  corresponding  nonspinning  case,  \diereas  the 

i 

[ initial  v-pro£iles  (Fig.  8)  were  obtained  by  superin^osition  of  the  non- 

spinning v-profiles  with  an  assumed  profile  to  account  for  the  spining 

Spinning  v = nonspinning  v + pure  spinning  v 


O<0<^ 


v<e<2n 


Fig.  8.  Initial  profile  for  p -direction  march. 

of  the  boundary.  The  assumed  profile  in  our  calciiLation  was  taken  to  be 

an  exponential  one.  With  initial  u-  and  v-profiles  so  obtained,  the 

corresponding  w-profiles  can  be  calculated.  Calculations  using  these 

initial  profiles  were  started  without  the  slightest  difficulty. 

rS  8 9 101 

In  the  pointed-body  problem \ initial  profiles  for  meridional 
march  near  the  tip  are  taken  to  be  the  nonspinning  ones  on  the  grounds 
that  since  at  the  tip,  the  body  radius  is  zero,  so  is  the  body's  spinning 
velocity.  The  fact  is,  of  course,  that  the  initial  station  must  always 
be  located  some  finite  distance  away  from  the  geometric  tip,  and  the  radius 
is  not  only  nonzero  but  also  may  even  be  not  small  as  in  a blunt  body  case. 

I 

I 

i 

I 

t 

L ______ 


F 


23 


Therefore,  although  this  method  is  reasonable  for  a pointed  body,  our 
method  with  a modification  for  pure  spinning  is  much  desirable  for  more 
general  bodies. 

3.4.b.  Initial  Profiles  for  Circumferential  Marching 

For  the  calculation  of  a new  crossplane  (fixed  in  y-plane) , scheme 
2 was  used  to  determine  the  first  station  (station  4 in  Fig.  9)  at  the 


Fig.  9.  Initial  profile  for  0 -direction  march. 

plane  0*0  using  the  solutions  at  known  stations  1-3,  5.  The  resulting 
solution  is  then  used  as  the  required  initial  profile  for  the  circum- 
ferential march.  Scheme  2 does  not  require  known  data  on  the  p-constant 
line  \diich  passes  through  station  4,  and  is  hence  especially  suitable 
for  the  present  purpose.  The  shaded  area  represents  the  dependence  zone 
of  station  4.  Evidently  the  mesh  of  scheme  2 does  not  con^jletely  enclose 
the  dependence  zone  uiless  Ap  is  taken  to  be  extremely  small.  The  latter 
is  certainly  inqjractical  because  the  con5)uting  time  will  become  exceed- 
ingly long.  However,  since  scheme  2 is  used  only  at  this  single  initial 
station,  incomplete  fulfillment  of  the  dependence  rule  at  a single  station 
is  inconsequential.  Some  effort  was  made  to  remedy  this,  but  it  was  found 
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later  unnecessary  and  hence  abandoned.  The  remedy  is  as  follows: 

As  the  calculation  was  carried  around  the  body  once,  the  solution  for 
9*0  was  compared  to  that  for  9 = 2n,  and  the  solution  of  9 * 2ir  was 
used  as  an  initial  profile  for  a second  round  of  calculations  between 
9*0  and  60®.  Only  small  differences  were  found  between  9*0  and 
9 * 27t  solutions  and  between  the  first  round  solution  and  the  second 
round  solution  for  0®<9l60®. 

3.5.  Ccmparison  with  Other  Methods 

In  the  scheme  used  by  Dwyer  and  Sanders  , the  crossflow  deriva- 
3 

tive  operator  (v  ^ in  our  notation)  is  always  kept  positive  by  eval- 
® 9 

uating  the  derivative  according  to  the  sign  of  the  velocity  v. 

Che  advantage  claimed  for  their  scheme  is  that  the  convection  process  is 
properly  modeled;  another  advantage  is  that  initial -values  are  not  needed 
in  the  same  fixed-u  plane  passing  the  unknown  station  being  calculated. 

Such  change  frcmi  forward  to  backward  differencing  according  to  the 
sign  of  V has  been  used  in  inviscid  flow,  and  indeed  models  the  convection 
properly.  What  should  be  differentiated  here  is  that  in  an  inviscid  flow, 
the  influence  and  dependence  are  determined  by  convection;  in  a three-dimen- 
sional boundary  layer  flow,  the  same  are  determined  by  diffusion  in  addition 
to  the  convection.  Diffusion  spreads  the  influence  or  dependence  across  the 
layer  instantaneously,  and  hence  smears  out  the  demarcation  between  the 
layers  of  opposing  flow-direction  from  idiat  would  occur  on  the  basis  of  con- 
vection alone.  The  upshot  is  that  this  technique  is  valid  in  the  inviscid 
case,  but  is  not  valid  vhen.  transplanted  to  the  boundary  layer  case. 
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Schematically  the  dependence  zone  for  the  station  AA2  CFig.  10)  is 

the  whole  wedge  AA2BB2DD2.  The  scheme  used  by  Dwyer  and  Sanders 

(23) 

ii?)lies  that  only  two  sub-wedges,  CCj^  and  A^A2C^C2B^B2  need 


Fig.  10.  Dependence  zone  for  Dwyer  and  Sanders'  scheme. 

be  accounted  for.  It  is  true  that  their  scheme  is  implicit  in  the 
normal  direction  along  which  all  grid  points  are  calculated  simultaneously; 
and  in  this  sense  the  diffusive  nature  is  loosely  modeled,  tfowever  it 
does  not  seem  that  an  implicit  calculation  would  svqjplant  the  require- 
ment of  satisfying  the  dependence  rule  cai?)letely. 

The  second  advantage  - initial  values  for  circumferential -direction 
march  are  not  needed  - would  be  a great  asset  for  a spinning -body  problem, 
but  such  advantage  is  seriously  offset  by  the  inability  to  satisfy  the 
dependence  rule  conpletely  unless  the  meridional  step-size  is  very  small. 

In  this  respect,  the  scheme  used  by  Dwyer  and  Sanders  has  the  identical 
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problan  as  our  scheme  2 discussed  before.  In  both  cases,  no  initial 
profiles  are  needed  in  the  lateral  direction  and,  as  shown  in  Fig.  9, 
there  vdll  always  be  an  area  vMch  cannot  be  covered  by  the  coiqputation 
mesh  unless  extrranely  small  steps  Au  are  taken.  The  difference  lies  in 
the  fact  that  in  our  case,  scheme  2 is  only  used  at  a single  point - 
station  at  9 = 0®,  while  Dwyer  and  Sanders  use  their  scheme  for  the 
entire  calculation.  These  authors  stated  that  the  rule  of  dependence 
zone  was  grossly  violated  in  their  calculation,  and  their  step-size  was 
10  times  larger  than  that  required.  They  further  remarked  that 
'•useful  results  were  obtained  \diich  agree  closely  with  the  converged 
solutions  that  obeyed  the  zone  of  influence.  One  possible  reason  for 
this  behavior  could  be  the  very  small  influence  of  spin  on  boundary  layer 
structure."  The  present  author  doubts  that  this  statement  is  generally 
valid. 

Kitchens  et  al^^^^  made  a con?jarative  study  of  four  different 
schemes,  including  our  earlier  version of  scheme  2 (Sect.  3.4b),  and 
the  scheme  of  Dwyer  and  Sanders;  calculations  were  made  for  a parabolic 
flow  over  a flat  plate.  Their  study  was  mainly  motivated  by  interest  in 
the  spinning-body  boxmdary  layer  problem.  The  scheme  they  proposed  in- 
corporated Lax's  difference  approximation  for  one -dimensional  unsteady 
flow.  This  seems  to  be  a good  feature  and  the  rule  of  dependence  zone  can 
be  completely  satisfied.  However,  as  with  Dwyer  and  Sanders'  scheme 
theirs  also  uses  one -step  in  the  u -direct ion,  and  tivo  steps  in  the  9- 
direction,  and,  again,  no  lateral  initial -profiles  are  required.  Our 
proceeding  comment  about  such  an  advantage  versus  the  inability  to  satisfy 
the  dependence  zone  for  the  Dwyer  and  Sanders'  scheme  hold  here  also.  In 


both  cases,  the  penalty  one  has  to  pay  far  outweighs  the  advantage 
received  in  return.  Secondly  Kitchens  et  al.  stress  the  importance 
of  their  scheme  being  more  insensitive  to  violations  of  the  zone  of 
dependence.  They  contend  that  such  insensitivity  "is  ia^xsrtant  be- 
cause.. .strict  adherence  to  the  zone  of  dependence  rule  could  be  very 
costly  because  of  the  small  downstream  step  sizes  required."  Actually 
the  necessity  of  using  small  downstream  step  sizes  arises  precisely 
from  the  adopted  scheme,  which  has  the  advantage  of  requiring  no  lateral 
initial  profiles.  In  our  approach,  such  lateral  initial -profiles 
are  generated,  as  computation  proceeds.  Only  at  the  initial  station 
(i.e.  9 * 0)  is  the  dependence  zone  (using  scheme  2)  not  con^jletely 
covered  by  the  con5)utation  mesh;  for  0 > 0,  the  Krause  schane  we  used 
could  readily  satisfy  the  dependence  zone  and  the  meridional  step  size 
is  not  subject  to  the  restriction  just  discussed.  Rather  than  stressing 
how  much  a particular  scheme  can  violate  the  dependence  zone,  it  seems 
preferable,  instead,  to  enphasize  how  such  violation  of  the  dependence 
rule  could  be  made  unnecessary. 

(91 

The  parabolicized  Navier- Stokes  approach  used  by  Lin  and  Rub  in  ^ 
also  does  not  require  initial  profiles  in  tte  circumferential -direction, 
but  this  is  no  bargain  because  treating  the  crossflow  in  an  elliptical 
fashion  is  a more  complicated  and  expensive  alternative. 
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4.  RESULTS  AND  DISCUSSION 

Two  cases  with  6"  and  30®  incidence  have  been  calculated  in 
this  work.  As  in  our  previous  non-spinning  problem^^^’^^^ , these  two 
cases  were  chosen  to  illustrate  problems  at  low  and  high  incidences. 

In  each  case,  results  for  two  different  spinning  rates  “ 0-5  and  1.0) 
were  obtained.  These  two  rates  correspond  to  the  maximum  spinning 
velocity  of  the  body  (v„)  being  less  and  greater  than  the  maximum 
circumferential  velocity  at  the  outer  edge  of  the  boundary  layer, 

For  example,  at  6®  incidence,  u = 0.5  and  1.0  gives  fVi,)  /V  =0.64 
and  1.28,  respectively,  (v^)^^  occurs  at  the  mid-body  where  the 
cross-sectional  diameter  is  maximum,  whereas  occurs  in  the  meridional 
planes  0 = 90®  and  270®.  The  results  of  6®  incidence  are  presented  in 
detail  in  Sect.  4.1  - 4.9;  the  results  of  30®  incidence  are  summarized 
in  Sect.  4.10. 

4.1.  Inviscid  Input 

The  inviscid  solution  for  the  present  spinning  case  is  the  same 
as  that  for  the  previous  nonspinning  case^^^'^^\  For  convenient 
reference,  the  distribution  of  sxorface  pressure  and  of  inviscid  velocities 
are  reproduced  in  Figs.  Il-rl2.  These  distributions  contain  the  synmetry  plane 
9 = 0 and  ir. 

4.2.  Boundary  Layer  Profiles 

The  roost  isqxjrtant  quantities  in  the  present  problem  are  the  cir- 
cumferential flow  profiles  which  take  on  a greater  variety  of  shapes  than 
is  usiially  encountered  in  any  single  problem  including  the  spinning  sharp 
cone  studied  before . Fig.  13a-e  illustrate  a few  examples.  Detailed 


plots  for  Figs.  13b-e  are  shown  in  Figs.  Bl(a)-(c)  of  Appendix  B. 
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Fig.  11. 


Fig.  12.  Surface  pressure,  a ■ 6*. 
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Fig.  13a  gives  the  v-pro£iles  .for  the  nonspinning  case'-  in- 
cluded here  for  the  purpose  o£  con^iarison.  A syinmetTy  plane  is  clearly 
indicated.  Reversed  flow  on  the  leeside  occurs  on  the  aft  body.  The 
pro£iles  are  non-dimensionalized  with  respect  to  local  outer-edge  velocity. 

Figs.  13b-d  show  the  v-pro£iles  with  the  spin  rate  oi  * O.S  at 
the  £ore-  (y  = - 0.80),  mid-  (y=  0)  and  a£t-body  (y  = 0.5).  For  the 
spinning  case,  it  is  more  convenient  to  non -dimens ionalize  the  pro£iles 
with  respect  to  Vg,  the  spinning  velocity  o£  the  body.  Vg  is  equal  to 
whg  and  hence  varies  longitudinally  as  h^  £or  a Fixed  03. 

Qualitatively  the  source  o£  these  proFile  shapes  can  be  visvialized 
by  super ijiqposing  the  nonspinning  proFile  on  a pure  spinning  proFile  (Fig.  8) 
and  by  comparing  the  velocities  Vg  and  V.  Consider  Fig.  13d,  For  example, 
the  proFile  For  0 = 20®  monotonically  decrease  From  the  body  toward  the 
outer  edge  with  Vg  being  greater  than  V.  That  For  9 = 90°  monotonically 
increases  with  V being  greater  than  Vg.  The  proFiles  oF  0 = 120°,  140° 
are  non-monotonic,  but  always  positive,  apparently  because  a strong  re- 
versed Flow  occurs  there  in  the  nonspinning  case  (Fig.  13a).  Such  reversed 
Flow  is  opposite  to  the  spinning  motion  and  hence  is  eFFectively  cancelled 
out.  The  proFiles  For0*  220°,  270°  are  non -monotonic  and  reversed  vdiile 
that  oF  0 » 340°  is  monotonic  and  reversed. 

In  view  oF  the  variety  oF  proFile  shapes,  we  have  to  use  the 
descriptive  words  '^nonotonic'*  and  ’’reversed"  to  distinguish  various  cases. 
"Reversed"  means  here  that  the  proFile  contains  both  positive  and  negative 
velocities.  In  an  upstream-moving  wall  problem  including  in  the 

present  problem,  a monotonic  reversed  proFile  is  Fairly  common,  but  a non- 
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mono tonic  reversed  one  Csudi  as  that  for  9 = 220®,  Fig.  13d)  is  not; 
the  latter  occurs  only  \diere  a strong  reversed  flow  already  exists, 
even  without  spinning.  Spinning  has  the  effect  of  reinforcing  such 
reversed  flow. 

Higher  spin  rates  exert  greater  influence  in  shaping  the  v- 
profile.  This  is  clearly  demonstrated  in  Fig.  13e  for  w = 1.0  at 
y “ 0.  Comparing  Fig.  13e  to  Fig.  13c,  it  is  evident  that  all  pro- 
files in  Fig.  13e  are  monotonically  decreasing  from  the  body  toward 
the  outer  edge,  whether  they  are  reversed  or  not. 

The  asymmetry  of  the  v-profiles  on  the  right  side  (0®  ' 180“) 
with  respect  to  those  on  the  left  side  (180“  - 360“)  is  the  major  source 
of  Magnus  effects.  The  difference  between  the  profile  shape  near  the 
body  on  two  sides  accounts  for  the  skin  frictional  contribution  to  the 
Magnus  force.  The  difference  between  the  overall  shape  of  those  pro- 
files on  two  sides  determines  the  centrifugal  pressure  contribution. 

The  meridional  velocity  profile,  u-profile,  does  not  show 
particularly  different  features  compared  to  the  corresponding  nonspinning 
ones,  although  a certain  degree  of  asymmetry  between  the  righv  and  left 
sides  is  also  noticeable  through  the  coupling  between  the  u-and  v- 
velocities.  A typical  san?)le  of  the  u-profiles  is  shown  in  Fig.  Bl(d). 

4.3.  Cross -Plane  Flow  Pattern 

The  cross-plane  flow  patterns  are  shown  in  Figs.  14a-e  for  co  ■ 1.0 
at  five  stations  extending  from  the  front  to  the  rear.  The  arrows 
indicate  the  flow  direction  determined  by  tan  ^ (hgW/v)  and  the  patterns 
are  presented  in  the  9, z -plane.  The  solid  curves  drawn  therefrom  represent 
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the  cross -plane  streamlines.  It  is  understood  that  these  streamlines 
are  not  the  actual  streamlines,  but  only  the  projection  of  three- 
dimensional  streamlines  on  the  cross -plane. 

The  most  obvious  effect  of  spinning  is  that  the  flow  near  the 
body  spins  with  the  body  due  to  viscosity,  so  that  the  streamlines 
(near  the  body)  are  parallel  or  nearly  parallel  to  the  body  surface. 

In  contrast,  the  outer  portion  of  the  boundary  layer  is  not  much  affected 
by  spinning  so  that  symmetry  with  respect  to  the  planes  0*0  and  it  is 
still  preserved. 

Ihe  structure  of  the  cross -plane  flow  can  be  better  analyzed  by 
studying  the  singularities.  These  are  located  at  the  intersection  of 
the  curves  along  i^ch  one  of  the  velocity  con^jonents,  v or  w becomes 
zero.  They  are  denoted  in  Figs.  14a  - e by  the  dotted  and  broken  lines. 
Before  we  discuss  Figs.  14a  - e further,  let's  start  with  the  non- 
spinning case  to  provide  a conparison. 

A typical  synnietrical  flow  pattern  for  the  nonspinning  case  is 
depicted  in  Fig.  15a.  The  solid  lines  with  arrows  are  the  streamlines 


Fig.  IS.  Crossplane  flow  pattern  (nonspinning) ; 
(a)  without  reversed  flow,  (b)  near  reversed 
region. 
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vMch  change  from  an  inflow  (toward  the  body)  to  an  outflow  (away  from 
the  body)  upon  crossing  the  line  w * 0.  Along  the  symmetry- lines  AA’, 

BB*,  CC*,  v-velocity  vanishes.  At  the  surface  A'C’,  both  v and  w 
vanish.  The  points  A*,  B’,  C are  the  saddle  points,  but  A',  C becomes 
a nodal  point  when  the  cross -plane  flow  becomes  overwhelmingly  an  outflow 
in  the  aft -body.  Also,  when  the  circxanferential  flow  reverses  on  the 
leeside,  a nodal  flow  pattern  develops  above  the  point  B*  (Fig.  ISb) . 

Now  we  return  to  the  spinning  case.  In  Fig.  14a  for  y = -0.90 
at  the  fore -body,  the  intersection  of  the  curves,  v = 0 and  w ■ 0, 
determines  a nodal  point  N and  a saddle  point  S.  The  nodal  pattern  near  N 
is  in  the  outer  portion  of  the  boundary  layer  and  is  only  slightly  affected 
by  spinning  so  that  N is  not  located  exactly  on  the  plane  0 = it  as  in  the 
nonspinning  case.  The  saddle  pattern  around  S is  located  in  the  middle 
of  the  fluid  flow  and  represents  a unique  feature  of  the  effect  of  spin. 

For  conqjarison,  one  may  think  of  the  curve  v * 0 in  Fig.  14a  corresponding 
to  the  line  BB’C’C  of  Fig.  15a  and  point  S of  Fig.  14a  corresponding  to 
point  B'  of  Fig.  15a.  In  other  words,  spinnning  motion  makes  the  flow 
asymmetric,  shifts  the  curve  v = 0 up,  and  moves  the  singularity  S away 
from  the  body  surface  into  the  fluid. 

At  y »-0.60  (Fig.  14b),  the  nodal  point  N disappears  and  the  curve 
w « 0 splits  into  two  branches.  On  the  other  hand,  the  cvirve  v = 0 and 
the  saddle  point  S are  displaced  upward  and  also  to  the  right.  Such 
trends  continue  up  to  y » 0.2  (Fig.  14c).  Here  the  two  branches  of  the 
curve  w * 0 are  further  separated,  vdiereas  the  curve  v * 0 and  the  point  S 
shift  farther  upward  to  the  right. 
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At  u » 0.4  in  the  aft -body  CFig*  14d) , the  tvro  branches  of  the 
curve  w * 0 join  together  at  the  left  upper  comer.  The  saddle  point 
disappears  altogether  and  the  flow  becomes  conqjletely  outward  Cw  > 0) 
except  at  the  left  upper  comer.  Tig.  14e  further  indicates  that  the 
curve  w = 0 conqiletely  disappears  at  y = 0.50  so  that  an  outflow 
prevails  over  the  entire  cross -plane.  Such  pure  outflow  would  have 
been  in^xjssible  in  a two-dimensional  cylinder  flow  from  the  viewpoint 
of  mass  continuity.  It  could  occur  here  because  of  the  accompanying 
meridional  (or  longitudinal)  flow,  a three-dimensional  effect. 

At  downstream  aft -body  stations,  the  effect  of  spinning  becomes 
more  and  more  pronounced.  ’:iiong  other  factors,  this  is  exhibited  by 
the  fact  that  the  layer  which  spins  with  the  body  becomes  thicker  and 
thicker  (Figs.  14d,  e).  The  layer  of  such  parallel  (or  nearly  parallel) 
flow  increases  especially  rapidly  between '0  = 180*  and  0 = 270*.  This 
is  because  the  flows  from  opposite  directions  collide  in  that  angular 
region.  In  particular,  the  reversed  flow  which  occurs  even  in  the  non- 
spinning case  is  reinforced  by  the  spinning  motion  there. 

The  present  problmn  constitutes  an  example  of  a three-dimensional 

moving -wall  boundary  layer  problem.  The  two-dimensional  moving  wall 

boundary  layer  problem  has  been  discussed  for  years  by  Sears  and  his 
f25") 

associates'-  ■'  who  classify  it  into  downstream-moving  (Fig.  16a) 
and  upstream-moving  (Fig.  16b)  cases.  Hiey  indicate  that  the  profile 
and  the  streamline  pattern  at  separation  for  the  upstream-moving  case 
might  appear  as  depicted  in  Figs.  16c,  d.  In  our  .present  problem,  the 
crossplane  flow  involves  both  a downstream-moving  (0-0-u)  and  an  upstream- 


Fig.  16.  Two-dimensional  moving-wall  problem; 
(a)  downstream-moving,  Cb)  upstream -moving , 
ic) , Cd)  profile  and  streamline  pattern  at 
separation  for  case  (b) . 


moving  (TT<0<2Tr)  wall  problem.  It  is  interestiiig  to  note  that  our 
calculation  indeed  deteimines  a saddle-point  pattern  for  Tt<0STr  in 
Figs.  14a-c,  although  not  much  can  be  said  about  the  connection  with 
respect  to  flow  separation.  The  separation  profile  of  Fig.  16c  has 
not  been  found  in  our  crossflow  results. 

The  above  crossplane  flow  patterns  may  be  redrawn  around  the 
body's  cross  sections  to  give  a more  physical  picture  (Fig.  17). 

Since  only  a limited  number  of  arrows  can  be  put  into  such  a pre- 
sentation, the  streamline  patterns  drawn  therefrom  are  necessarily  less 
accurate,  llie  previoas  nonspinning  patterns  are  also  included  for  com- 
parison. Since  the  latter  is  syiimetrical,  only  half  of  the  cross  section 


is  shown. 


The  results  presented  thus  far  are  for  u » 1.0.  For  the  lower  spin 


rate  (u  ■ O.S),  the  general  trends  are  similar  and  the  layer  which  spins 
with  the  body  is  less  noticeable.  The  pattern  for  (u*1.0aty  = 0.2 
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looks  similar  to  that  for  u = 0.5  at  p = 0.4,  implying  a lag  of 
the  spinning  effect.  The  patterns  for  u » 0.5  are  shown  in  Figs.  B2 
and  B3. 

4.4.  Meridional -plane  Flow  Pattern 

In  conjunction  with  the  preceding  crossplane  patterns,  the  flow 
patterns  in  fixed  meridional  planes  (0  ■ 0®,  90®,  180®,  and  270®)  are 

also  plotted  in  Figs.  18a  - d.  The  flow  direction  is  given  by 
tan'^(h^w/u)  and  is  represented  by  the  arrows  in  y-z  coordinates.  . A 
curve  along  \diich  w ■ 0 could  have  been  added  C^ut  not  done  here)  in 
Figs.  18a -d  to  separate  area  of  inflow  (with  respect  to  the  body)  from 
that  of  outflow.  The  calculated  longitudinal  flow  is  always  pointed 
downstream  (i.e.  u >0).  From  the  fore-body  toward  the  aft-body,  the 
flow  generally  turns  increasingly  outward.  The  meridional -plane  patterns 
shown  are  less  con^jlicated  than  the  preceding  crossplane  patterns.  Decrease 
of  u)  tends  to  alleviate  the  outflow  in  the  aft -body,  but  otherwise  similar 
patterns  prevail.  The  details  are  shown  in  Figs.  B4.  The  meridional - 
and  cross -plane  flow  patterns  together  help  one  to  visualize  the  three- 
dimensional  flow  structure,  although  it  is  difficult  to  depict  the  stream- 
line patterns  in  three -dimensions.  IVhen  both  the  cross-  and  meridional- 
plane  flows  become  strong  outflows,  idien  the  resultant  flow  also  must 
be  a strong  outflow  leading  inevitably  to  the  beginning  of  flow  separation. 

4.5  Meridional  (Or  Longitudinal)  Skin  Friction,  c^^ 

The  meridional  skin  friction  distribution  is  shown  in  Fig.  19a 
for  01  ■ 0.5.  This  distribution  is  syiiinetric  over  the  fore  body,  but  be- 
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comes  gradually  asymmetric  further  aft.  A conqjarison  for  c^^  at 
the  mid-body  Cu  “ 0)  for  different  spin  rates  is  shown  in  Fig.  19b. 

Higher  spin  rate  makes  tlie  distribution  even  more  asymmetric.  How- 
ever, the  contribution  from  the  asynmetry  of  c^  tp  the  Magnxis  force 

is  negligible. 

The  minimura-c^  curve  shown  in  Fig.  19a  indicates  that  this 
ccmqjonent  of  skin  friction  becomes  vanishingly  small  near  u = 0.55  and 
0 ■ 310®.  This  marks  the  beginning  of  the  reversal  of  the  meridional 
flow.  Since  the  circumferential  flow  is  already  reversed,  our  calcu- 

. i 

lation  must  be  terminated  idienever  meridional  reversal  starts.  Other 
common  symptoms  associated  with  flow  separation  were  also  noted  at  the 
juncture,  so  that  the  reversal  of  meridional  flow  is  taken  here  to  be 
the  onset  of  separation.  What  happens  is  that  starting  from  9 ■ 310® 
at  u ” 0.55,  the  separation  spreads  quickly  in  the  circumferential 
direction  as  indicated  by  the  broken  line.  The  distribution  of  c^^  for 
0)  ■ 1.0  shows  little  difference  and  is  given  in  Fig.  B5. 

4. "6  Circumferaitial  Cor  Crossflow)  Skin  Friction,  c^g 

The  variation  of  c^g  over  the  body  for  tu  ■ 0.5  is  shown  in  Fig.  20a  I 

which  can  be  better  understood  by  referring  to  the  v-profiles  in  Figs.  13b -d. 

Over  the  fore-body  (y  ■ -0.90),  c^g  is  positive  for  9 <180®  (except  near 
9-0®),  but  negative  for  0 >180®.  At  downstream  stations  (say  y * -0.60), 

I 

c^g  begins  to  be  negative  for  9 < 130® . Farther  downstream  (y  > 0) , C£g  | 

becomes  positive  for  180®  < 0 < 360®.  This  change  of  sign  can  be  explained 

by  con5)aring  the  profiles  for  9 - 210®  in  Figs.  13c, d.  In  one  case  (Fig.  13c), 
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a monotoixic  reversed  profile  results  in  a negative  c^g,  in  the  other 
CFig.  13d) , a non-monotonic  reversed  profile  gives  a positive  c^g  . 

The  reason  for  the  latter  profile  was  discussed  in  Sec.  4.2.  The  large 
positive  C£g  in  the  area  just  discussed  ijnplies  strong  reversed  cross- 
flows and  is  related  to  the  eventual  flow- separation.  Con5)lete  distri- 
butions of  C£y  and  C£g  over  a spinning  body  (such  as  Figs.  19a  and  20a) 

have  not  been  available  before.  For  the  sharp  cone  problem,  calculated 

f8  9") 

distributions  of  skin  frictions  are  either  completely  absent^  ’ ^ or 
given  only  for  a single  crossplane. 

Fig.  20b  shows  the  effect  of  spin  rate  on  C£g  along  the  mid- 
section, u=  0.  The  increase  of  asymmetry  due  to  higher  spin  rate  is 
clearly  indicated.  The  distribution  for  u = 0.5  falls  roughly  midway 
between  those  for  0 and  1.0.  Such  an  almost  linear  relationship 
apparently  reflects  the  fact  that  the  circumferential  flow  near  the  body  is 
controlled  by  the  spinning  motion.  Fig.  B-6  gives  the  variation  of 
over  the  body  foru*  1.0.. 

4.7.  Centrifugal  Pressure,  p^^ 

The  circumferential  flow  always  generates  a centrifugal  pressure,  ! 

but  this  induced  pressure  causes  no  concern  in  the  nonspinning  case  because  | 

i 

the  contributions  from  the  two  sides  of  the  symmetry-plane  cancel  each  | 

other.  The  asyranetry  of  the  circumferential  flow  profiles  CSec.  4.2)  due 
to  the  body’s  spinning  leads  to  an  asymmetrical  distribution  of  the  centri- 
fugal pressure,  p^^  as  shown  in  Fig.  21a  for  u * 0.5.  In  the  circumferential 
direction,  the  variation  of  Ip^^l  follows  roughly  that  of  the  inviscid  v-velocity 
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Fijz.  21.  Centrifugal  pressure,  (a)  Distribution  of 
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(Fig.  11);  maximum  )p^^|  occurs  at  9 * 90®  where  the  v-velocity 
across  the  boundary  layer  reaches  a maximum.  In  the  meridional  direction, 
increases  toward  the  rear  as  the  boundary  layer  thickens.  Ip^^l  in- 
creases with  the  spin  rate  as  is  evident  in  Fig.  21b,  idiich  shows  a com- 
parison of  p^^  along  the  mid-section,  y = 0,  between  oj  = 0.5  and  oj  = 1.0. 

A similar  general  trend  of  p^^-distribution  for  u = 1.0  is  given  in  Fig.  B7. 

4.8.  Displacement  Thicknesses 

4.8a.  v-Velocity  Displacement  Thickness,  A* 

9 

Spinning  introduced  certain  new  features  in  Ag  . First  it  is  noted 
that  A?  can  become  negative.  This  arises  because  the  flow  velocity  on  the 
body,  Vg,  can  be  greater  than  the  outer  velocity  V (Fig.  22b).  Thus  the 
defining  integral 


9 Jo 


(V-v)dz 


becomes  negative.  Physically  this  means  that  viscosity  causes  a mass  flow 
excess  rather  than  defect,  and  the  effective  body  boundary  is  displaced 
into  the  body  rather  than  into  the  fluid.  However,  the  idea  of  negative 
displacement  thickness  is  only  of  conceptual  significance  here.  Since  the 
circumferential -flow  displacement  thickness  is  imch  smaller  than  that  of 
the  meridional -flow,  the  total  displacement  thickness  is  still  positive. 
Secondly  the  above  integral  becomes  indeterminate  when  V vanishes. 

The  latter  happens  at  9 » 0 and  ir.  In  the  nonspinning  case,  vanishing  V 

* 

leads  to  the  vanishing  of  the  v-profile  across  the  boundary  layer,  so  Ag 
can  be  simply  set  to  zero.  In  the  spinning  case,  tMs  is  not  true,  q is 
chosen  to  replace  V as  the  normalizing  factor  (see  Eqs.  9a,b).  The  latter 
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j * *1 

( introdxxes  a complication  in  that  Aq  becomes  negative  for  n<Q<2v  ,1 

I " 

■ vrfiereV  is  negative.  The  negative  A.  discussed  above  for  positive  V 

i ^ 

\ has  physical  meaning;  the  negative  A here  does  not.  This  awkward 

situation  coiHd  have  been  avoided  if  one  added  a minus  sign  in  front  of 
; the  integral  in  Eq.  9c;  but  then  A^  would  become  discontinuous  at  9 * 180®. 

Fig.  22a  gives  the  distribution  of  A*  for  u » 0.5.  The  maximum 
magnitude  in  the  9-direction  occurs  around  9 = 270®.  This  is  to  be  expected 
physically  because  of  strong  reversed  flow  there.  Fig.  22c  coiqjares  A* 
for  various  spin  rates  at  the  mid-section,  u = 0,  with  the  higher  spin  rate 
(ta  » 1.0)  causing  Ag  to  be  negative  throughout.  To  be  consistent,  A*  for 
the  nonspinning  case  is  considered  in  Fig.  22c  to  be  negative  for  ir£B-2Tr 

although  conventionally  this  is  taken  to  be  positive.  A similar  distribution 

* 

of  Ag  for  10  “ 1 ; 0 is  given  in  Fig . B8 . 

* 

4.8b.  Total  Displacement  Thickness,  A 

The  total  displacement  A*  is  calculated  using  Eq.  (n) . Its  distribution 

$ 

for  (0  “ 0.5  is  shown  in  Fig.  23a.  At  the  present  low  incidence  (a  » 6®)  , 
the  meridional  flow  dominates  over  the  crossflow,  and  the  total  displace- 
ment thickness  A*differs  very  little  from  the  corresponding  u-velocity  dis- 
placement thickness  (Fig.  23b).  At  the  fore-body.  A*  is  small  and  nearly 
uniform  around  the  body.  Farther  aft.  A*  increases  and  the  maximum  thickness 
position  shifts  from  180®  into  the  direction  of  spinning,  so  that  the  dis- 
placed body  is  twisted  longitudinally.  Fig.  23a  shows  slight  irregularities 
near  9*0®  and  9 » 180®  over  the  aft-body.  This  is  a computational  effect, 
but  does  not  affect  the  rest  of  the  results.  The  asymraetr)'  of  A*  di«tribution 
again  reflects  the  spinning  effect  on  the  longitudinal  flow  through  its 
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coupling  with  the  crossflow.  A similar  distribution  of  A*  and  A* 
for  w = 1.0  are  shown  in  Figs.  B9,  BIO. 

The  momentum  thicknesses  defined  by  Eq.lOa-d  were  also  calcu- 
lated but  not  presented  here. 

4.9.  Magnus  Force 

Majoi  contributions  to  the  Magnus  force  come  from  (1)  circum- 
ferential skin  friction,  (2)  centrifugal  pressure,  and  (3)  the  asymmetric 
displacement  effect.  The  last  factor  has  not  been  determined  here  be- 
cause a separate,  inviscid-flow  program  is  needed  to  calculate  the 
presstire  over  the  displaced  body.  The  magnitude  of  the  other  two  and  the 
resultant  Magnus  force  are  calculated  lising  Eqs.  12b,  13c,  and  14;  the  results 

I 

are  shown  in  Figs.  24a,  b for  oj  = 0.5  and  1.0,  respectively.  The  frictional 
and  the  centrifugal  side-forces  c^^  and  c^^,  partly  compensate  each  other. 

At  lower  spin  rate,  w = 0.5,  the  frictional  force,  c^^,  is  larger  than 
the  centrifugal  force,  c^^  (Fig.  24a),  resulting  in  a positive  Magnus  force, 

Cmg  to  th'e  right  in  Fig.  24a).  At  higher  spin  rate,  u = 1.0, 

the  opposite  is  obtained — i.e.,  the  centrifugal  force  is  larger  in  magni- 
tude than  the  frictional  force,  and  the  Magnus  force  becomes  negative  (i.e., 
to  the  left  in  Fig.  24b) . 

The  change  of  the  Magnus  force  from  positive  at  lower  spin  rate 
to  negative  at  higher  spin  rate  has  been  previously  observed  experimentally^^^^ . 
This  change  was  ascribed  there  to  wake  vortices.  Incidentally,  Ref.  26  pro- 
vides an  extensive  literature  of  experiments  on  the  Magnus  effect  not 
mentioned  in  the  present  work.  A series  of  flow-visualization  pictures 
can  be  found  in  Ref.  27. 


Fig.  24.  Sectional  side -force. 
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4.10.  The  High- Incidence  Problem 

In  our  previous  nonspinning  high- incidence  work^^^’^^^,  calcu- 
lations in  the  circumferential  direction  were  carried  out  only  to 
the  open  separation  lines.  Determination  of  the  leeside  separated  flow 
(shaded  in  Fig.  25a)  was  attenqjted  later,  but  the  results  have  never 
been  reported.  We  have  included  that  effort  here  as  a part  of  the 
extension  to  the  spinning,  high- incidence  problan. 

For  the  nonspinning  case,  calculation  was  extended  into  the 
region  where  the  crossflow  reverses.  The  resulting  distributions  for 
the  skin  frictions  are  shown  in  Figs.  25b, c.  The  meridional  skin  friction 
C£^  , shows  a sharp  decrease  upon  entering  the  reversed  region  from  the 
wliidward  direction  and  then  increases  somewhat  toward  the  leeside 
symmetry -plane  0 = 180®  (see  p = -0.88  to  u = -0.86).  Such  a trend  in- 
tensifies so  rapidly  that  c^.^  becomes  zero  at  y = -0.833  and  negative 
thereafter  near  6 =»  140®.  The  calculation  was  therefore  terminated. 

To  put  the  matter  into  the  right  prospective,  we  should  recall 

fl2'l 

that  in  the  nonspinning  case  for  6®  incidence,'-  C£^  was  similarly 
found  to  decrease  sanehwat  faster  upon  entering  the  reversed  flow  region 
from  the  windward  side,  but  the  variation  in  the  longitudinal  direction 
was  so  gradual  that  C£^  did  not  become  zero  except  on  the  aft -body  and 
the  reversed  flow  region  extending  over  a quarter  of  the  body  length  was 
calculated. 

The  quick  termination  at  high  incidence  leads  us  to  conclude  that 
the  leeside  separated  region  associated  with  an  open-separation  can  not 
be  determined  according  to  the  classical  boundary  layer  theory.  This  is 
by  no  means  an  obvious  conclusion  from  the  outset,  because  it  does  not 


seem  to  be  certain  that  the  flow  involved  violates  in  any  way  the 
boundary  layer  assumptions,  at  least  not  in  the  iqjstream  fore-body. 

Our  experience  in  the  low  incidence  case^^^^  just  mentioned  especially 
encourages  us  to  think  this  way.  In  any  case,  the  early  termination 
for  the  high  incidence  is  somewhat  unexpected,  and  it  is  unclear  vdiat 
the  precise  reasons  are  for  it  although  one  can  always  cite  that  the 
reversed  flow  is  too  strong  and  the  boundary  layer  assim5)tions  such 
as  the  potential  pressure  are  no  longer  valid. 

Our  interests  in  extending  the  high  incidence  problem  to  include 
the  spinning  morion  are  multifold:  Ci)  to  determine  vrtiether  boundary 

layer  calculation  can  be  started  even  though  there  is  no  rigorous  way 

» 

to  obtain  the  initial  profiles,  (ii)  to  find  out  whether  spinning  motion 
would  change  the  crossflow  sufficiently  so  as  to  be  able  to  extend  the 
calculation  further  downstream,  and  Ciii)  to  examine  the  flow  features 
from  solutions  which  can  be  obtained  by  the  present  suiqjler  method,  even 
for  a limited  region  of  the  body  surface. 

In  our  previous  non-spinning  high- incidence  work, the 
problem  involved  first  the  calculation  of  stagnation  region,  then  the 
symmetry -plane  solution,  the  nose  region  solution,  and  finally  the  main- 
body  solution.  Spinning  couples  these  local  flows  together  and  hence 
they  cannot  be  determined  separately  in  the  order  described.  Mb  rational 
way  is  now  known  to  resolve  this  question. 

In  the  nonspinning  case,  different  starting  methods  were  used  for 
the  low-and  high- incidence  problems.  For  the  high- incidence  case,  the 
starting  method  eii5)loyed  the  nose  region  solution  as  just  mentioned.  For 


the  low  incidence  case,  approximate  profiles  were  made  up  on  the 
jxistification  that  the  boundary  layer  varies  only  slightly  from  the 
windside  toward  the  leeside  over  the  fore-body. 


For  the  spinning  case,  the  starting  profiles  for  low  incidence 
(6®)  were  described  in  Sect.  3.4a.  The  question  then  is  whether  the 
same  technique  can  be  enployed  for  the  high- incidence  case.  Mien  this 
was  tried  out,  it  was  found  (someidiat  surprisingly)  that  calculation 
proceeded  downstream  without  difficulty.  The  starting  station  was 
located  at  y = -0.943.  This  is  not  for  any  other  reason  than  that  the 
nonspinning  solution  for  this  particular  y- station  previously  calculated 
was  still  available  on  tape.  This  nonspinning  solution  modified  by  a 
pure  spinning  profile  (Fig.  8)  makes  up  the  required  initial  profiles. 

Two  spinning  cases  with  u ■ 0.5  and  1.0  are  calculated  for  the  same 
ellipsoid  (Fig.  2)  at  30®  incidence  Sample  results  of  the  v-velocity  pro- 
file are  shown  in  Fig.  26  for  oi  ■ 0.5  at  y = -0.836.  More  details  of  these 
profiles  are  given  in  Fig.  Bll. 
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■Die  calculated  meridional  skin  friction,  c^,  for  o ■ SO® 
and  (0  “ 0.5  is  shown  in  Fig.  27a.  Con^jared  to  Fig.  19  for  6*  incidence, 
the  c^ -distribution  along  a fixed  U”Station  changes  noticeably  within 
a short  longitudinal  distance  (fran  y ■ -0.940  to  y « -0.833  in  Fig.  27a). 
Along  y =■  -0.87  or  -0.85,  C£^  increases  from  9 ■ 0*  to  6 =•  90®^/  100®, 
decreases  sharply  around  0 = 140®  160®,  then  increases  again  slightly 

tp  to  9 » 220®  '\i240®.  Around  0 ® 250®,  a second  dip  appears.  At  y = 
-0.833,  c^^  drops  almost  to  zero  at  9 = 160®  so  that  the  present  calcu- 
lation must  be  terminated.  This  occurs  at  the  same  y-station  as  for 
the  nonspinning  case.  The  corresponding  distribution  of  c^g  is  shown 
in  Fig.  27b,  vMch  displays  similar  trends  as  Fig.  20a. 

The  skin  frictions  for  co  « 1.0,  a = 30®  are  given  in  Fig.  B-12. 
Increase  of  the  spin  rate  enables  the  calculation  to  be  extended  downstream 
slightly  (terminated  at  y » -0.812),  but  otherwise  the  same  characteristics 
prevail. 

Thus  the  question  of  initial  profiles  did  not  present  an  unsurmount - 
able  problem,  but  calculation  could  not  be  continued  in  the  longitudinal 
direction  very  far  after  the  point  at  which  an  open  separation  for  the 
corresponding  nonspinning  case  occurs.  Since  the  starting  point  of  open 
separation  continues  to  move  upstream  as  the  incidence  increases'*  this 
means  that  the  area  over  vhich  the  spinning  boundary  layer  can  be  calcu- 
lated decreases  at  higher  incidences. 
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5.  CONCLUSIONS 

The  three-dimensional  laminar  boundary  layer  over  a spinning 
blunt  body  of  revolution  at  incidence  is  investigated  here. 

Complete,  concrete  results  are  obtained  for  an  ellipsoid  of  revolution 
vdiich  provides  a model  of  more  general  geometry  in  comparison  with  the 
sharp  cone  studied  before.  Spinning  introduces  unique  features  from 
the  viewpoints  of  both  fluid  mechanics  and  methodology. 

Spinning  causes  asymmetry  which,  in  turn,  results  in  the  side 
forces  known  as  Magnus  effects.  The  asynmetry  is  most  pronounced  in 
the  crossflow,  but  the  meridional  flow  is  also  affected  through  the 
coupling  as  evidenced  in  the  asynmetrical  distribution  of  c^^  and 
Ilie  variety  of  crossflow  profiles  found  is  greater  than  usually  en- 
countered in  any  single  flow  problan.  The  calculated  flow  patterns 
in  the  cross -and  meridional -planes  provide  a clear  picture  of  the  flow 
structure.  In  particular,  the  cross-plane  flow  patterns  and  the  singu- 
larities involved  are  presented  in  detail,  along  with  the  skin  frictions 
and  displacement  thicknesses.  A new  derivation  of  the  total  displacement 
thickness  is  given  and  a negative  displacement  thickness  is  found  to  be 
physically  meaningful  in  the  spinning-body  boundary  layer,  although  its 
effect  on  the  overall  flow  is  negligible  in  the  case  calculated.  The 
Magnus  forces  due  to  the  circumferential  skin  friction  and  the  centri- 
fugal pressure  are  determined;  these  two  forces  partly  condensate  each 
other.  At  lower  spin  rates,  the  frictional  force  is  larger,  resulting 
in  a positive  Magnus  force.  At  higher  spin  rates,  the  opposite  is  the 
case.  The  overall  magnitude  of  the  Magnus  force  increases  with  spin  rate. 
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At  high  incidence  C30°),  boundary  layer  calculations  are  carried 
out  in  the  longitudinal  direction  only  up  to  vdiere  an  open  separation 
occurs.  The  leeside  separated  region  associated  with  an  open  separ- 
ation is  not  amenable  to  a classical  boundary  layer  treatment.  It  is 
not  certain  whether  this  difficulty  can  be  overcome  merely  by  using  an 
experimental  pressure  distribution  instead  of  the  potential -flow  pressure 
used  here. 

Methodologically,  the  present  work  strictly  follows  the  boundary 
layer  theory,  vduch  does  not  become  invalid  merely  because  of  the 
presence  of  the  crossflow  reversal.  The  associated  zone  of  dependence 
can  be  effectively  observed  with  the  Krause  scheme.  The  generation  of 
initial  profiles  to  start  the  calculation  did  not  lead  to  an  impasse 
as  previously  conceived.  Profiles  for  the  longitudinal  march  can  be 
approximated  by  stqjerin^xjsing  the  corresponding  nonspinning  profile 
with  an  assumed  pure  spinning  profile.  Those  for  the  circumferential 
march  are  generated  as  the  con5)utation  proceeds.  The  present  calcu- 
lation terminates  idien  the  meridional  skin  friction,  c^^  , becomes  zero 
because  the  boundary  layer  theory  ceases  to  be  valid  when  both  the  u- 
and  the  v-profiles  reverse.  In  this  respect,  the  parabolicized  approach 
offers  no  advantage.  To  save  conpiting  time,  experiments  were  made  to 
show  that  larger  step-sizes  can  be  used,  depending  on  how  severe  the 
relevant  adverse  pressure  gradient  is.  The  case  of  u * 1.0  at  b"  inci- 
dence was  calculated  with  less  than  one  hour  of  CPU  time  on  a UNIVAC  1108 
conpater. 
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APPENDIX  A 


Derivation  of  the  Equation  Determining  the  Total 
Displacanent  Thickness  for  Three-Dimensional 
Boundary  Layer 
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In  analogy  vdth  the  usual  reasoning  in  deriving  the  mass  continuity 
equation,  we  consider 
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Net  change  of  the  mass  defect  across  the  side-surfaces 


3^  ^ 

3x  3y 


(A-2) 


which  is,  of  course,  balanced  by  the  mass  flow  across  the  top  surface 

of  the  volume  element,  but  this  is  not  of  our  concern  here. 

Next  we  consider  the  mass  defect  of  the  same  volume  element  in 

terras  of  a single  displacement  thickness  6 '’Fig.  AlCb))  which  has  no  simple 

integral  definition  as  Eq.  (A-1)  for  6 and  5 . 6 is  so  defined  that  the 

X y 

net  mass  defect  across  the  side -surfaces  for  the  same  volume  element  in 

terras  of  5 will  be  the  same  as  Eq.  (A-2)  in  terms  6 and  5 . The  latter 

X y 
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condition  is  required  in  order  to  ensure  the  mass  continuity  upon  replacing 

5 and  5 by  5,  while  the  mass  flow  condition  at  the  top  surface  is  kept 
X y 

unchanged. 

Following  the  same  steps  by  comparing  the  outgoing  and  the  incoming 
mass  defects,  we  obtain  the  net  change  of  mass  defect  across  the  side-surfaces 
in  terms  of  5 

_ ^ ^ ^ _ (A3) 

3x  3y 

Equating  Eq.  (A- 3)  to  Eq.  CA-2)  yields 


3U6xhy  ^ ^ 5yhx 
3x  * 3y 


3U6hy  3V6hx 
3x  * 3y 


9 


or 

> |y[VC6-V]’  CA4) 

which  is  the  equation  first  derived  by  Moore. To  specialize  this 
equation  to  the  streamline  coordinates,  we  may  set 


Then  Eq.  (A-4)  becomes 
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Integrating 

hyU(6-5^) 


X 

0 


+ 


|^Ch^U6j,)dx 


0. 


(A*S) 


Assume  that  x = 0 is  the  stagnation  point  where  the  streamline 
originates  andU  = 0.  Also  in  a streamline  coordinate,  ® , hence 

Eq.(A-5)  becomes 

^ ^ 

vdiich  is  the  form  given  by  Lighthill . Lighthill  derived  Eq.  (A-6)  by 
four  different  methods,  and  refers  to  Moore's  derivation  as  the  method  of 
velocity  con^arison.  Our  derivation  follows  the  reasoning  employed  for  the 
mass  continuity  equation,  and  is  someidiat  akin  to  Lighthill 's  flow  reduction 
method,  but  differs  in  the  formalism.  Lighthill 's  derivation  is  also 
specialized  to  the  streamline  coordinates.  It  is  straightforward  to  include 
the  compressibility  effect  in  Eq.  (A-4)  or  to  write  it  in  slight  by  different 
form  as  Eq.  (11)  of  the  text. 
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APPENDIX  B 
Additional  Figures 


This  appendix  contains  a number  of  additional  figures  obtained  from 
this  investigation. 


Fig. 

Bl. 

Fig. 

B2. 

Fig. 

B3. 

Fig. 

B4. 

Fig- 

B5. 

Fig. 

B6. 

Fig. 

B7. 

Fig. 

B8. 

Fig. 

B9. 

Fig. 

BIO 

Fig. 

Bll 

Fig. 

B12 

v-velocity  profiles,  a = 6®,  u = 0.5  and  1.0. 

Cross -plane  flow  pattern  (rectangular),  a * 6®, 
u » 0.5  and  1.0. 

Cross -plane  flow  pattern  (circular),  a » 6®, 

(D  » 0.5  and  1.0. 

Meridional -plane  flow  pattern,  a=  6®,  oj  = 0.5. 

Meridional  skin  friction,  a = 6®,  w * 1.0. 

Circumferential  skin  friction,  a * 6®,  u = 1.0. 

Centrifugal  pressure,  a = 6®,  u = 1.0. 

v-velocity  displacement  thickness,  a = 6®,  u = 1.0. 

u-velocity  displacsnent  thickness,  a = 6®,  u * 1.0. 

Total  displacement  thickness,  a * 6®,  u = 1.0. 

v-velocity  profiles,  a » 3.0®,  u * 0.5. 

Skin  frictions  for  a = 30®,  u = 1.0 

Ca)  Cb) 


i 


Fig.  Bl(d).  u-velcx:ity  profile. 


0.0482  y = -0.80  Az  = 0.0482  y = -0.60 
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Fig.  B2.  Crossplane  flow  pattern  (rectangular) 
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Fig.  B2.  Continued 
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Fig.  B2.  Crossplane  flow  pattern  (rectangular) 
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Fig.  B3.  Crossplane  flow  pattern  (circular),  a » 6® 
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